Dynamics of Distillability 
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The time evolution of a maximally entangled bipartite systems is presented in this paper. The 
distillability criterion is given in terms of Kraus operators. Using the criterion, we discuss the 
distillability of 2 x 2 and n x n(n > 2) systems in their evolution process. There are two distinguished 
processes, dissipation and decoherence, which may destroy the distillability. We discuss the effects 
of those processes on distillability in details. 
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Maximally entangled states are an essential ingredi- 
ent in most applications of quantum information[l,2]. In 
particular, in quantum communications, by using entan- 
gled states several proposals were devised to transmit 
secret messages between two locally separated parties [3- 
5]. These proposals have been successfully implemented 
experimentally by using pairs of photons generated via 
parametric down conversion[6-9]. On the other hand, 
quantum entanglement is a fragile feature, which can 
be destroyed by interaction with the environment. To 
overcome this difficulty, Bennett etal., Deutsch etal. and 
Gisin etal. presented several schemes to distill maximally 
entangled states of two qubits out of a set of pairs in 
mixed entangled states[10-12]. These process are called 
entanglement distillation, which only consists of local ac- 
tions and classical communication. Unfortunately, it is 
not known in general, which mixed state are distillable at 
that time. Later on, the Horodecki family shown that all 
entangled states of 2 x 2 and 2x3 systems can be distilled 
into a singlet [13], and they proved that nonpositivity of 
the partial transposition is a necessary condition for the 
state of an arbitrary bipartite system to be distillable[14]. 
But, it does not be the sufficient one. In fact, there are 
states which have non-positive partial transpose [15] but 
they are not distillable [16- 18]. Recently, a useful suffi- 
cient criterion, the so-called reduction criterion, has been 
derived[18]. It shown that one can explicitly construct a 
protocol to distill p if there exists some vector satis- 
fying 



(V>|Tr B p®l-/#) <0, 



(1) 



where Trg stand for the trace over the second subsystem. 
Moreover, some results of the distillability may be gen- 
eralized to the case of continuous variable systems [19]. 

Since the early study on the entanglement and distill- 
ability, most works in this area are mainly concentrated 
on separability and distillability for a concrete class of 
mixed states. In practice, however, the destruction of a 
maximally entangled state is closely related to a dynam- 
ical process. One of the examples is that the destruction 
of a maximally entangled state due to the interaction 



with the environment is a dynamical process. Some au- 
thors become aware of the importance of the dynamical 
properties for quantum entanglement [20, 21]. 

In this paper we investigate the time evolution of the 
distillability for a initially maximally entangled state. We 
consider a simple bipartite system which consists of two 
particles with the same dimensions. Two cases are taken 
into account in this paper, the first case is that only 
one of the two particles is subjected to the environment, 
the other case is two particles in the bipartite system 
are entirely under the effects of the environment. From 
the viewpoint of dimensions, this paper may be divided 
into two parts, i.e., 2x2 systems and n x n (n > 2) 
systems. For the 2x2 systems, there is an alternative 
necessary and sufficient condition for distillability. For 
high dimension, however, we only have necessary or suf- 
ficient criterion for distillability, but not a criterion for 
both. Our results show that the decoherence (caused by 
the environment considered here) do not change the dis- 
tillability of the 2x2 systems, whereas the distillability 
for a dissipation systems remains unchanged only within 
a short time scales. For high dimension, however, the 
distillability condition is more complicated. It does not 
only depend on the coupling of the systems to the envi- 
ronment, but also on the initial condition. 

Distillability. -We consider two parties, Alice and Bob, 
who share several pairs of particles. For simplicity, we 
assume here that the particle has the same dimension d. 
Each pair is initially in a maximally entangled state 



1 d 



(2) 



Interactions with the environment transform this pure 
state into mixed state. This process may be described 
by a linear, trace-preserving, completely positive map L 
as 



L (Po) = Pf =^2Aip A\, 



(3) 



where po stands for the initial state, throughout this 
paper we assume that po — \ip + )(^j + \,i.e., the system 



1 



is initially in a maximally entangled state. The trace- 
preserving property implies that the operator A4 obey 
the constraint 



(4) 



with 1 the identity matrix on the Hilbert space. One of 
the physical implementations of the Kraus operators Ai 
is as follows. We consider a system interacting with the 
environment. The evolution of the total system (system 
plus the environment) is govern by a unitary operator 
U(t). The reduced density matrix of the system may 
be given by tracing the total density operator over the 
environment, 



p f (t) = Tr B [U(t) Po (0) ® \0 B ) (Ob | U f (t)] , 



(5) 



where |0s) stands for the initial state of the environment. 
For a set of complete and orthonormal bases {p. B } of the 
environment, one obtains 



Pf(t) = J2 Ai P° A i> 



(6) 



with Ai = {/j, B \U(t)\0B} is the so-called Kraus operator. 
It is evident that Ai satisfy Y^Jl A\Ai = 1) f° r U (t) is uni- 
tary. So L(p ) — > pf(t) is a trace preserving completely 
positive map. 

For 2x2 systems, the sufficient and necessary condition 
for distillability is 



Ff = Tr(p p f ) > -, (7) 
In terms of the Kraus operators, this condition is 

u 

1 



> 



(8) 



In derivation of eq.(8), we used the initial condition 
po = \i/j + )(ip + \. For high dimension (> 2), there is a re- 
duction criterion, it states that if there exists some vector 
lib) such that 



Gf = {ip\Tr B Pf ® 1 - p\i/>) < 0, 



(9) 



then the final state pf is distillable. An important aspect 
of this criterion is that if one finds a state \ip) satisfying 
eq.(9), then one can explicitly construct a protocol to 
distill p/[18]. For a initially maximally entangled state 
, we suppose that the most promising state, by which 
the distillation proposal is constructed, is \ip + ) itself. In 
practice, we may always perform a distillation before the 
maximally entangled state decoheres far away from its 
initial state. In this sense the sufficient condition for dis- 
tillability is 



G f = £<V> + |Tr s (A i |V> + )<V + l4) ® !|^ + > 



i=l 

k 

<^KV> + MMV> + >| 2 . 

i=l 



(10) 



Therefore if we get the Kraus operators, we may know 
exactly the distillability of a state. Some words of cau- 
tion arc now in order. The Kraus operators arc not 
unique in general. For example, in eq.(6) we may choose 
the other bases {v B } instead of {p B } to compute the 
Kraus operators. However, different sets of Kraus opera- 
tors, which describe the same dynamical process, may be 
transformed each to other by a unitary transformation. 
In this sense, the distillability criterion do not depend on 
the choice of Kraus operators. 

2x2 system with decoherence.- Decoherence occurs 
due to unwanted interactions between our quantum sys- 
tem and its environment. These interactions cause only 
information leak out of the system. Typically, this pro- 
cess may be described by the following two Kraus oper- 
ators [22]. 





-ft I > 



A, 
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VT - e 



-27* 



(11) 



We consider a pair of entangled particles a and b. If only 
one of them (say a) are subjected to environment, the 
time evolution of a maximally entangled state 



\lP + ) = ^(\0aA h }-\l ai h )) 



is then given by 



Pf(t) = 






Vo 
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(12) 



where 7 a is the decay rate for particle a. This process 
may occur when a photon from a entangled pair is trans- 
mitted through a fiber whose length is randomly mod- 
ulated by acoustic waves, or an atom from the pair is 
exposed to interactions with a environment that consists 
of a set of harmonic oscillators [24]. Substituting eq.(ll) 
into eq.(8), we obtain 



Ff = 



1 1 



-7a t 



(13) 



It is evident that the final state is always distillable. If 
the two particles are entirely under the effect of the envi- 
ronment, it is easy to show that the new Kraus operator 
for the whole system are 



Ai = Ai a ® A u , A 2 = A 



2a 



■A 



2b: 



(14) 
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So, Ft in this case is 

F f = \ + \e-^^\ (15) 

the final state is always distillable, too. 

2x2 system with dissipation. - Different from the case 
of decoherence, dissipation leads not only to the decay of 
off-diagonal elements of the density matrix, but also the 
energy loss (decay of the diagonal clement of the density 
matrix). The effect of energy loss to the environment 
is usually described by a master equation[25]. which in 
the Born-Markov approximation may be represented in 
terms of Kraus operators, 

*-(; MS 

If we only transmit one of the entangled particles through 
a noisy channel, the Ff is then given by 

2 

F /= ^|(V; + |A|V + )| 2 = e-^. (17) 

z=l 

It is greater than \ only for t < 2ln2/j a . This indi- 
cates that we must perform distillation within 2ln2/j a 
in order to distill a maximally entangled particles. Sim- 
ilarly, we obtain Ff = e s -6- for the case that the 
two entangled particles are both under the effects of the 
environment. It is well known that the character time of 
decoherence is much shorter than the dissipation one, so 
we can ignore the effects of dissipation in general. From 
the viewpoint of distillation, however, the dissipation is 
more destructive, for the information loss due to deco- 
herence may be reconstructed by distillation, but the in- 
formation loss caused by dissipation does not. Figure 
1 shows a results of Ff(t) defined in eq.(7) by solving 
master equations 

p = -i[H , p] + - ^ faiPPi - oj<Jlp - po+crr), 

i—a^b 

and 

p = -i[H ,p] + ( 2a iP a i ~ a t a iP - P a i a i) 

i—a : b 

numerically, where Hq represents the free Hamiltonian 
and is defined as Hq = ui(a* + a^). It is well known that 
the first master equation describes the dynamics for a 
dissipation system, while the second one governs a deco- 
herence process. Both master equations can be derived 
by using Markov approximation and assuming system- 
bath interaction J2i( a+a i + h.c.) and ^(o^cij + h.c), 
respectively. Here we use a notation a; to denote the 
bath mode annihilation operator. It is clear that Ff(t) 
corresponding a dissipation process goes down below 0.5 



earlier than a decoherence process. In this sense, we say 
dissipation is more destructive. 
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FIG. 1. Ff(t) versus t, the dashed line is drawn for a deco- 
herence system, while the solid line is for a dissipation system. 
Parameters chosen are 7 — 0.6u>, t is chosen in units of 1/ui. 

n x n system with decoherence.- Without loss of gen- 
erality, we present the following model to describe the 
effect of decoherence, the model Hamiltonian is 

H = H s + Hb + Hi , 

Hi = J2fj({os})®9j({oB}), (18) 

3 

where fj and gj are functions of system operators and 
bath operators, respectively. In order to ensure that only 
the decoherence occur in the system, we make a con- 
straint [o s ,H s ] = on the operator o s . This condition 
implies that we may factorize the time evolution operator 
of the composite system in the following way 

N 

U(t) = e -KH,+H B )t . -Q e F 3 (t)-f 3 ({o s })-G 3 ({o B }) ^ ( 19 ) 
j=l 

Here. 

/i({o.}) • G 3 -({o fl })(3 = 1,2, ...,#)} 

are elements of Lie algebra enlarged by H s , Hb, Hi , while 
function Fj(t) is determined by 

?L e FAt>fj({°s}>Gi({oB}) = 

dt 

= F J (i)F 1 J (t)/ J ({ 0s })G({ 0s })e^( t )^« =»^« -». (20) 

Eq.(20) indicates that in bases spanned by the eigenstates 
of H s , the diagonal elements of the system density matrix 
remain unchanged in the evolution process, while the off- 
diagonal parts gain a complex phase with the evolution. 
A simple calculation give the element of reduced density 
matrix for the system 



3 



(m\p s \n) = Tr B (m\p\n) 

= Tr B {m\U(t)p s {0) $ ^(O)^ (t)|n) 

JV 

= TrB[e~ i(Emt ~ Ent) Y[ e Fl{t)IimGl({oB]) 

3=1 

• p m „(0)®p B (0)e F ^" G ^»] 
= p ro „(0)e-T-W- ir -«, (21) 

where / Jm satisfies fj{{o s })\m) = fj m \ m )- The quanti- 
ties 7 mn (t) have the following properties, 

7mn(t) = 0,for m = n, and 7 m „(t) > 

for otherwise. Physically, no energy transfer between 
the system and the environment require "f mn — for 
to = n, while 7 mn (t) > for m ^ n. For r mn (t), we 
have r mn (i) = —T* nm {t), this property is directly from 
the hermity of pf. For a initially maximally entangled 
state 



= ^^2\m a ,m b )(n a ,n b \ 



It is easy to check that 
G f (t) = Tr(Tr bPf ® 1 - p f ) = - £ -« 

• cos(r^(t) + r^(t)). 



(22) 



-7^„(*)-7m„(*) 



(23) 



Here, we assume the two entangled particles are both 
under the effect of the environment, 7^ n (t) and r^, n (t) 
(7mn(*) an( i rtmW) are defined by eq.(21) for a- 
particle(b-particle). Eq.(23) is a damping-oscillation 
function of time. If cos(r^ n (t) + T^ nn (t)) > 0, the final 
state of pf is distillable. If pf does not violate the Peres 
separability criterion, the final state pf can not be dis- 
tilled, this can be done for a given state p/. For a model 
presented in [23], r^ n (i), (x = a,b) may be rewritten 
explicitly as 



Q 2 



(m 2 - n 2 )(ut - sinujt)p(uj)duj, (24) 



where g u>x is the coupling constant of particle x to the 
environment, p(uS) stands for the spectrum distribution 
of the environment. Eq.(23) shows that the final state pf 
is not always distillable, its distillability would depend on 
the detailed information of the system even in the case 
of decoherence. 

n x n system with dissipation.- We consider a spin-n 
(d = 2n + 1) particle interacting with its environment. 
Under the Born-Markov approximation, the dissipation 
process may be described by a master equation in the 
Lindblad form[25] 

p = -i[H ,p] + l(2S_pS + - S + S- P - pS+S-), (25) 



where S + (S-) is spin operator defined by S+ = (S x +iS y ) 
(S- = (S x — iS y )). Hq = £IS Z is the free Hamiltonian 
of the spin-n particle. 7 is the decay rate. We would 
like to note that eq.(25) only described the time evolu- 
tion of one particle in the entangled pair. In general, the 
particles in the pair experience different environment, for 
they are transmitted trough different noisy channels. By 
using the method presented in[26], we obtain the element 
of the density matrix up to the first order of 7 

1 jt 



Pm„(t) = ^ + -!jy/(f + to +!)(/- to)(/ + n + !)(/ - n) 



it—(ilm + fin — i~f{f + to)(/ — to + 1) 
ij(f + n)(f- n+1)). 



(26) 



In derivation of eq.(26), the initial condition p(0) = 
2 Em n l TO )( n l was used- If the two particles are both 
exposed to the interactions of the environment, follow- 
ing the procedure presented in case A of this section, we 
arrive at 



1 d , 

G /W = E 2Re ^ + lty/U + m + l)(f-m) 



■ y/(f + n + l)(f-n) 
— it(Qm + Qn — vy(f + ni)(f — m + 1) 
-il(f + n){f-n+l))] 2 . 



(27) 



According to the sufficient condition(8), the final state is 
distillable if G/(t) < 0. For 7 = 0L2f2, the dependence 
of the critical time t c defined by Gf(t c ) = on dimen- 
sion of the systme is shown in Fig. 2. As Figure 2 shows, 
the larger the dimension of the system, the shorter the 
character time of distillability destruction. 



FIG. 2. The dependance of the critical time t c on system 
dimension 
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Figure 3 is plotted for a numerical simulation of Gf (t) 
for a decoherence system and a dissipation system, the 
dynamics of the dissipation system is govern by master 
equation (25), while the time evolution of the decoher- 
ence system are described by 

p = -i[H , p] + ~(2S z pS z - S z S z p - pS z S z ), (28) 

0.1 I 1 1 1 1 1 1 1 , 1 



o - ( 




Time t 



FIG. 3. Gf(t) as a function of time. Solid line is for a deco- 
hence system with dimension 5, dotted line for a dissipation 
one with dimension 7, dashed line for decoherence with di- 
mension 7 and dotted and dashed line for dissipation with 
dimension 5. In this plot, we choose 7 = 0.6S1 and time t is 
plotted in units of 1/Q 

From figure 3 we can see that in the case of dissipation 
Gf(t) go up to zero after a short time evolution from our 
initial state IV'(O)) = 4 J2 m l m ' — m ); whereas Gf(t) for a 
decoherence process is always below zero. We would like 
to address that this conclusion is not a general one, which 
would depend on the initial state of the system and the 
bath information, for example, if we choose m = 0, and 
p(uj) — constant in eq.(24), we may find a time when 
cos(T a mn (t)+T b mn (t)) <0. _ 

In summary, the dynamics of distillability for a bi- 
partite system are investigated in this paper. The de- 
struction of the maximal entanglement for the system is 
closely related to the interactions with the environment. 
There are two kinds of interactions which lead to the de- 
struction of maximal entanglement. One is the quantum 
decoherence, the another is the dissipation. From the 
viewpoint of distillability, the dissipation is more harmful 
in quantum communication. For example, the decoher- 
ence do not change the distillability of a 2 x 2 systems, 
whereas the dissipation does. For high dimension, both 
dissipation and decoherence do destroy the quantum en- 
tanglement, but in different ways. 

We would like to note that the initial state consid- 
ered here is a maximally entangled state. If the initial 
state belongs to a special class of entangled states, the 
local environment can enhance the quantum entangle- 



ment [27] from the view point of quantum teleportation. 
This increases the efficiency of distillation and makes the 
undistillable state to be distillable . In the framework 
of quantum information theory, the state change allowed 
by quantum mechanics may be classified into three types. 
The first one is the unitary evolution, the second is the 
interaction with a environment, and the last one is a mea- 
surement performed on the quantum system. The uni- 
tary evolution is of cause change the state of the bipartite 
system, for the two particles in the bipartite systems may 
interact each other. Although we do not discuss the dy- 
namic of distillability of such a system in this paper, the 
method of this paper can easily generalized to this case. 
As to the measurement, the method presented here is 
also available, because the most general type of measure- 
ment can be understood within the framework of unitary 
evolution[28]. In fact most generalized measurement can 
be realized through many dynamical processes [29]. 
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